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The generating function for recurrent Polya walks on the four 
dimensional hypercubic lattice is expressed as a Kampe-de-Feriet 
function. Various properties of the associated walks are enumer- 
ated. 

Lattice statistics play an important role in many areas of chemistry and 
statistical physics. Random lattice walks have been investigated extensively, 
both numerically and analytically, and explicit analytic expressions for recur- 
rent lattice walks (those that end where they began) are available for several 
lattices in one, two and three dimensions [1,2]. The interest in this prob- 
lem follows from the numerous connections between random walks and other 
problems in mathematical physics. They form the basis of certain proofs in 
the problem of self-avoiding walks [3], they are a key component in the study 
of lattice dynamics [4], and they have an extensive following in the mathe- 
matical literature [1,3]. Furthermore, d-dimensional lattice Green functions 
have recently [5] been shown to underlie the theory of d-dimensional stair- 
case polygons and to be related to the generating function for d-dimensional 
multinomial coefficients. 

The aim of this note is to extend this list to four dimensions by summa- 
rizing a study of the lattice Green function 

P(*) = P&) = ~J I [ r ; TTTy \A < 1 (1) 
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which generates the statistics for recurrent Polya walks ( uniform step prob- 
abilities) on the 4d hypercubic lattice. 

From (1) it is clear that P(z) is an analytic function of z 2 within the 
unit circle and has a branch point at z 2 = 1. It is representable as a single 
integral in a variety of ways. From the relations in [6], for example, we have 



and 



where 



P(z) = J e-'li(-xz)dx (2) 
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kl = -[l± x 2 z 2 ] Jl - -x 2 z 2 - (1 - -x 2 z 2 )Vl-x 2 z 2 ] (4) 

and K denotes the complete elliptic integral of the first kind. 

In terms of the variables u = z 2 and p(u) = P{z), P(z) satisfies the 
Fuchsian differential equation 

4w 3 (w-4)(w-1)p ( ™ ) (w)+8w 2 (5m 2 -20u+l2)p'"(u)+u(99u 2 -293u+U2)p"(u) 

+ (57u 2 - 106m + 16)p'(u) + (3u - 2)p(u) = (5) 
Asymptotic analysis of (2) and series analysis shows that near z 2 = 1 

P(z) = Co + d(l- z 2 )\n(l -z 2 ) 

C = 1.23946712... (6) 

d = 2/tt 2 

Upon expanding (1) in powers of z 2 , one has 

oo 

P(Z) = E *nZ 2n (7) 
n=0 



where 
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On the other hand, by expanding the Bessel function in(2), integrating term- 
by-term, and proceeding as in the appendix in [6], we find 

n r jVgW Fl -n/2, (l-n)/2, -n, 1/2; 

a " = ^H 4Fs[ 1/2 -n, 1/2 -n, 1; 1J ' (9) 

Next, by using the identity [6] 

4 F 3 (-n, (1 - n)/2, -n/2; 1/2 - n, 1/2 - n, 1; 1) = 

4 F 3 (-n, 1/2, 1/2, 1/2; 1/2 - n, 1/2 - n, 1/2 - n; 1) (10) 
expressing the latter 4 F 3 in terms of its series, and noting that 

(-n) fc = (-l) fc n!/(n - fc)!, (1/2 - n) fc = (-l) fc (l/2) n /(l/2) n _ fe , (11) 

we arrive at the double hypergeometric series 

- (l/2) m+TO (l/2) T 3 t (l/2)^ ^ 2 " 
() " m to (l)S^n n!m! (12) 

which defines a Kampe-de Feriet function [7]: 

P(z) = F£**[ 1,1,1; -; -; z\z 2 } (13) 

=« : 2 1 [ 1/ !'. 1; l 7 ?: ^74, .74]. 



The quantity a n in (8) is the probability that the random walker ends at 
its starting point in 2n steps (not necessarily for the first time). The number 
of such walks is therefore r n = 8 2n a n , which obeys the recursion relation 

n 4 r n - 4(20n 4 - 40n 3 + 33n 2 - 13n + 2)r n _ x 

+256(4n 4 - 16n 3 + 23n 2 - 14n + 3)r„_ 2 = 0. (14) 

The first ten values are given in the table. 

Number of 2n-step recurrent walks 
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n 


r„ 


1 


8 


2 


168 


3 


5120 


4 


190120 


5 


7939008 


6 


357713664 


7 


16993726464 


8 


839358285480 


9 


42714450658880 



10 2225741588095168 
For large n, for the d-dimensional hypercubic lattice [8] 

rW ~ 2(2d) 2n (d/Ann) d/2 

so in the present case a n ~ 2/n 2 n 2 and r n ~ 64 n (2/n 2 7r 2 ). (For example, this 
yields a 50 o ~ 8.10569 x 10~ 7 , compared to the exact value a 50 o = 8.0976 x 
10~ 7 , a difference of less than 0.1%). 

To our knowledge, the only previous studies of P±{z) were by Bender et 
al.[8], who examined PaiX) as a function of d and by a \jd expansion obtained 
a value consistent with (6), and by Mayer and Nuttall [9], who derived a 
differential equation consistent with (5). We have not examined P(z) outside 
the unit circle, although it is clearly analytic in the z 2 plane cut along the real 
axis from z 2 = 1 to infinity. While it is doubtful that P(z) can be expressed 
in terms of elliptic integrals as for the three dimensional cubic lattices [2], 
it is hoped that further analysis of (13) will lead to a more transparent 
representation of its analytic properties. 
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